Role of the two- and three-body hyperon-nucleon interaction in A-hypernuclei 
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We discuss the effect of the hyperon-nucieon-nucleon interaction on the binding energy of a A 
hyperon in closed shell hypernuclei from A=5 to 91. The A binding energy has been calculated 
using the Auxiliary Field Diffusion Monte Carlo method that we have used, for the first time, to 
study light and heavy hypernuclei within the same model. Our results show that including a three- 
body component in the hyperon-nucleon interaction leads to a saturation of the A binding energy 
remarkably close to experimental data. In contrast, the two-body force alone gives a very unphysical 
binding energy. 



PACS numbers: 21.80.+a, 26.60.Kp, 21.60.De 

The onset of strange baryons (E~ and A) in neu- 
tron matter at densities of order (2 — 3)/9o, where po — 
0.16 fm~ 3 , has been questioned for a long time. Recent 
theoretical calculations based on the Brueckncr-Hartree- 
Fock theory suggest that any process generating new on- 
shell degrees of freedom in high density Fermionic matter 
leads to a substantial softening of its equation of state 
(EoS) (see for example [1, 2] and references therein). 
When occurring in the inner core of a neutron star, such a 
mechanism would reduce the value of its predicted max- 
imum mass and of its radius. Until a few years ago, 
astrophysical observations of neutron stars were concen- 
trated in a relatively narrow region in the neighborhood 
of the Chandrasekhar limit (M ~ 1.41 M ). Most of the 
realistic EoS based on the hypothesis of matter made of 
nucleons only, while compatible with these observations, 
predict a maximum mass typically larger than 2 Mq. 
This result can be considered very robust. As an exam- 
ple, recent Quantum Monte Carlo (QMC) calculations of 
the equation of state of Pure Neutron Matter (PNM), 
Symmetric Nuclear Matter (SNM) and baryonic matter 
at /3- and /i-equilibrium using realistic Density Depen- 
dent Potentials (DDP) [3], essentially confirm the behav- 
ior predicted by Akmal, Pandharipande and Ravenhall 
with a full AV18+three-body interaction [4]. With such 
nuclear Hamiltonian the predicted EoS supports a max- 
imum neutron star mass larger than 1.97 Mq recently 
observed [5]. 

In so far as the appearance of strange baryons is con- 
cerned, the situation is more controversial. Some au- 
thors (see e.g. [6, 7]), suggest that the appearance of hy- 
perons in the EoS does not lead to very strong effects. 
Other recent papers, like Refs. [1, 2, 8-10], show a more 
substantial influence, but with contradictory outcomes in 
terms of the predicted maximum mass of neutron stars 
not compatible with the observations [11]. Therefore, the 
situation is far from being completely settled. 

A combination of reasons leads to the uncertainty in 
the analysis of the influence of strangeness degrees of free- 
dom in the EoS. First of all, the interaction between nu- 



cleons and hyperons is still far from being known with 
sufficient accuracy. The prospective measurements of 
properties of light hypernuclei should improve the qual- 
ity of the available data, making possible a realistic phc- 
nomenological analysis. Second, the theoretical tools em- 
ployed are all affected by uncontrollable intrinsic approx- 
imations as soon as one tries to push the study beyond 
few body systems. As a consequence, so far it is not 
clear how well the model hyperon-nucleon (YN) poten- 
tials work in the limit of medium mass hypernuclei, and, 
as a consequence, in the extrapolation to homogeneous 
matter. However, in the last few years important ad- 
vances have been made both on the experimental and on 
the theoretical side. Several experiments aim to measure 
the binding energy of different A-hypernuclei [12-14]. On 
the theoretical side, the development of Quantum Monte 
Carlo methods opened the way to study consistently nu- 
clear systems from few-nucleons to infinite matter [15-17] 
within the same scheme/model. 

In this letter we discuss the use the Auxiliary Field Dif- 
fusion Monte Carlo ( AFDMC) , to study a non-relativistic 
Hamiltonian based on a phenomenological AN interac- 
tion in order to show how the inclusion of explicit ANN 
terms provides the necessary repulsion to realistically de- 
scribe the separation energy of a A hyperon in hypernu- 
clei of intermediate masses. This point makes very clear 
the fact that the lack of an accurate Hamiltonian might 
be responsible for the unrealistic predictions of the EoS 
that would tend to rule out the appearance of strange 
baryons in high density matter. 

Several models can describe the hyperon-nucleon (YN) 
interaction. A number of potentials in the Nijmegen soft- 
core form have been developed in the past (like NSC89 
and NSC97x). A recent review of these interactions, to- 
gether with HF calculations have been published by Dapo 
et al. [18]. These potentials are accurate in describing the 
available scattering two body data, and have been used 
in the BHF studies previously quoted. Starting in the 
1980s, a class of Argonne-like interactions have been de- 
veloped by Bodmer, Usmani and Carlson on the grounds 
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of Quantum Monte Carlo calculations. A possible form 
of a three-body YNN interaction was also proposed in the 
same context [19-23]. More recently Polinder et al. [24] 
proposed a potential based on a chiral perturbation the- 
ory expansion. As an alternative a cluster model to study 
light hypernuclei has been recently proposed by Hiyama 
and collaborators (see for example Refs. [25, 26]). Inter- 
esting results on A-hypernuclei have also been obtained 
within a A-nucleus potential model, in which the need of 
a functional with a more than linear density dependence 
was shown, suggesting the importance of a many-body 
interaction [27] . 

An important point that needs to be made is that AN 
and ANN interactions are both at the two-pion exchange 
order. Another important difference with respect the nu- 
clconic case is that the mass of the intermediate excited 
state £ compared to the A is much smaller than in the 
pure nucleonic case, where the difference between the nu- 
cleon and the A is much larger. AN and ANN interactions 
should therefore be considered necessary in any consis- 
tent theoretical calculation. In 2002 Nogga et al. [28] per- 
formed Fadeev-Yakubowsky calculations of the + and 
the 1 + state of A H and A He in order to study charge 
symmetry breaking effects. In both cases they predict a 
A-separation energy that is too small, and claim that an 
attractive ANN interaction is necessary. 

We have revisited the problem from a slightly different 
starting point. We employed a potential in coordinates 
space, including an explicit repulsion between baryons, 
with NN, AN and ANN components. Keeping the pa- 
rameters of the AN interaction at the values determined 
by Usmani et al., we computed the ground state energy 
of a set of hypernuclei, and calculate for each the quan- 
tity B A , i.e. the separation energy of the A hyperon, by 
means of the AFDMC method, using a realistic nucleon- 
nucleon interaction. We select one of the possible set of 
parameters of the ANN interaction suggested in Bodmer 
et al., and then by Usmani and collaborators, that rea- 
sonably reproduces experimental data on a set of light 
hypernuclei. 

Within this model, nuclei and hypernuclei are de- 
scribed as non-relativistic particles interacting via two- 
and three-body forces: 
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Here A refers to the total number of baryons, nucleons 
plus the A particle. In order to test the effect of using dif- 
ferent nuclear Hamiltonians on the A-separation energy, 
we use three different two-nucleon potentials : the Ar- 
gonne AV4' and AV6' [29], that are simplified versions 
of the Argonne AV18 [30] and the Minnesota potential 



from Ref. [31]. 

Isospin conservation implies that a A hyperon can ex- 
change a pion only with a A7r£ vertex. This fact leads 
to the consequence that one pion exchange (OPE) pro- 
cesses are not allowed. The lowest order AN coupling 
must therefore involve the exchange of two pions, with 
the formation of a virtual £ hyperon, as illustrated in 
Fig. 1(a). One pion exchange processes can only occur 
through the exchange of a K, K* kaon pair. This pro- 
cess has the effect of exchanging the strangeness between 
the two baryons, as shown in Fig. 1(b). The AN interac- 
tion can therefore be modeled with a central term, which 
includes the AN exchange operator e(P x — 1), plus a spin- 
dependent contribution: 

VAi = v Q (r A i)(l-e + £P x ) + -v„T^(m. K r A i)cr A -cri, (3) 

where P x , vo and T% are defined in Ref. [23] and refer- 
ences therein, and a A and er^ are Pauli matrices acting 
on the A and nucleons. Both the spin-dependent and the 
central terms contain the usual tensor operator T v act- 
ing twice. All the pion exchange interaction is therefore 
active at intermediate range. The short range contribu- 
tions are as usual included by means of a phenomenolog- 
ical central repulsive factor, included in vo(r). For more 
details see for example Ref. [23]. 

The remaining diagrams in Fig. 1, are two-nucleon one- 
hyperon interactions, which are at the same TPE order, 
and should therefore be included together with the two- 
body part in order to have a consistent description. The 
three-body potential v A ij can be conveniently decom- 
posed in a contribution that we label as = v Ai A+v Ai p 
and that corresponds to the p-wave and s-wave two-pion 
exchange diagrams (respectively 1(c) and 1(d) in Fig. 1), 
and a dispersive term that includes short range contribu- 
tions, labeled as v Ai j. They can be expressed as: 
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The definition of the functions Xi A and Z(x) as well as 
the range of parameters for the three-body force can be 
found in [23] and references therein. 

The ground-state energy of the many-body nuclear and 
hypernuclear Hamiltonians, is computed by means of the 
AFDMC method. The algorithm was originally intro- 
duced by Schmidt and Fantoni [32] in order to deal in an 
efficient way with spin-dependent Hamiltonians. A trial 
wave function is propagated in imaginary-time r by 
sampling configurations of the system in coordinate-spin- 
isospin space. Expectation values are computed averag- 
ing over the sampled configurations. In the r — > oo limit, 
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Figure 1. Meson exchange processes between nucleons and hyperons. 1(a) and 1(b) represent the AN channels. 1(c), 1(d) and 
1(e) are the three-body ANN channels included in the potential by Usmani et al. [21, 23]. 



the evolved state approaches the ground-state of H and 
thus the ground state properties of the system can be 
obtained. 

For a system with A nucleons, the quadratic operato- 
rial structure of the nuclear Hamiltonians leads to 
a number of spin-isospin states in the propagated wave 
function which grows exponentially with A. This num- 
ber quickly becomes intractable as A gets large. Stan- 
dard Green's Function Monte Carlo (GFMC) calcula- 
tions are in fact limited to up to 12 nucleons [33] or 16 
neutrons [17]. By applying the Hubbard-Stratonovich 
transformation the computational cost of the calculation 
becomes proportional to A 3 and systems with a larger 
number of particles can be studied [16]. The AFDMC 
algorithm can be applied to nuclear systems interacting 
via the Argonne V6-type potentials, for which the two- 
body force can be separated into a spin-independent and 
a spin-dependent part. The latter can be written as a 
sum of real matrices which contain proper combinations 
of the components of V6. By means of the diagonaliza- 
tion of such matrices it is possible to write the imaginary 
time propagator in the Hubbard-Stratonovich form (see 
Refs. [16, 34, 35] for a detailed discussion). 

A straightforward variant of AFDMC can be applied to 
A-hypernuclear systems, including the two-body (3) and 
three-body (4) hyperon-nucleon interactions. It is indeed 
possible to recast the AN and ANN interactions so that 
they contain at most two-body operators. These terms 
can directly be included in the AFDMC propagator. The 
rest of the algorithm closely follows the nucleon-only ver- 
sion [16]. 

We assume that the wave function of a single A- 
hypernucleus is a nuclear Slater determinant (the same 
of Ref . [35] ) , multiplied by a single particle wave function 
for the A hyperon. As nucleon single particle states we 
use the radial solutions of the Hartree Fock problem with 
the Skyrme force and we consider a lsi/2 single particle 
state for the A particle. With the wave function defined 
we consider nucleons and the hyperon as distinct parti- 
cles. In this way, we do not include the AN exchange 
term of the AN potential directly in the AFDMC prop- 



agator, because it mixes hyperon and nucleon states. A 
perturbative treatment of this factor is however possible. 

A direct comparison of energy calculations with ex- 
perimental results is given for the A-separation energy, 
defined as: 

Ba = B nuc — Bhyp, (5) 

where B nuc and Bh yp are respectively the total binding 
energies of a nucleus with A nucleons and the correspond- 
ing hypernucleus with A nucleons plus one A. The most 
significant outcome of the calculation is the fact that the 
inclusion of the three-body ANN interaction qualitatively 
changes the saturation properties of the A-separation en- 
ergy. However, this result might depend on the particu- 
lar choice of the NN interaction used to describe both the 
nucleus and the hypernucleus. In particular, one might 
expect a strong influence from the different nucleon den- 
sity generated by disparate models. In order to discuss 
this possible dependence, we have performed calculations 
with different NN interactions having very different sat- 
uration properties. The nuclear Hamiltonians considered 
here are semi-realistic, and can be easily implemented 
within the AFDMC scheme. We should point out that 
in neither case did we use a three nucleon interaction. 

In Tab. I we show the results of the AFDMC sim- 
ulations for the A-separation energy in |He and ^O. 
For each hypernucleus, the two columns correspond to 
calculations using the AN interaction only or both the 
AN+ANN force of Ref. [23] with different NN interac- 
tions. As it can be seen, for ^He the extrapolated val- 
ues of Ba with the two-body AN interaction alone are 
about 10% off, and well outside statistical errors. In con- 
trast the inclusion of the three-body ANN force gives a 
similar A binding energy independently to the choice of 
the NN force. On the grounds of this observation, we 
feel confident that the use of AV4', that makes AFDMC 
calculations less expensive and more stable, will in any 
case return realistic estimates of Ba for larger masses 
when including the ANN interaction. We have checked 
this assumption performing simulations in ^O, where 
the discrepancy between the A-separation energy com- 
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puted using the different NN interactions and the full 
AN+ANN force is less than few per cent (last column 
of Tab. I). The various NN forces considered here are 
quite different. The AV6' includes a tensor force, while 
AV4' and Minnesota have a simpler structure. We com- 
pared the AV4' and Minnesota, which have a similar op- 
eratorial structure but very different intermediate- and 
short-range correlations. The fact that the inclusion of 
the ANN force does not depend too much on the nuclear 
Hamiltonian is quite remarkable, because the different 
NN forces produce a quite different saturation point for 
the nuclear matter EoS, suggesting that our results are 
pretty robust. 
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Table I. A-separation energies (in MeV) for ^He and ob- 
tained using different nucleon potentials (AV4', AV6', Min- 
nesota) and different hyperon-nucleon interaction (two-body 
alone and two- plus three-body). In the last line the experi- 
mental values of Ba are from Ref. [36] and Ref. [37]. 



The results on the A separation energies are summa- 
rized in Fig. 2. We compare the prediction of the hyperon 
binding energy in the AV4'+AN and AV4'+AN+ANN 
models for a few closed shell hypernuclei with the exper- 
imental values observed in the same mass range. While 
the results for lighter hypernuclei might be inconclusive 
in terms of the physical consistency of the ANN contri- 
bution to the hyperon binding energy, the computations 
for A \Ga and 9 ^Zr reveal a completely different picture. 
The saturation binding energy provided by the AN force 
alone is completely unrealistic, while the inclusion of the 
ANN force gives results that are qualitatively much closer 
to the experimental behavior. We should notice that 
the results might be further improved by a refitting of 
the terms in the ANN force. In particular, according 
to Ref. [21], in the present calculations the s-wave con- 
tribution is not present. Moreover, we are missing the 
explicit inclusion of the kaon exchange term. This con- 
tribution (see Eq. (3)) can be estimated at first order in 
perturbation theory by computing the expectation of the 
corresponding term. As an example, the values of the cor- 
rection on the total energy we obtained for e = 0.1 [23] is 
-0.33(6) MeV in A He and +0.2(4) MeV in ^O, the latter 
negligible compared to the corresponding binding energy. 

For 9 A Zr we should also consider a charge symmetry 
breaking (CSB) potential. The latter can be easily in- 



cluded as a term of the form 

v%* = T?CZ aB TZ(m„r Ai ), (6) 

amounting to an isospin dependent correction to the cen- 
tral potential. The inclusion of the CSB term using per- 
turbation theory would be zero in isospin-symmetric nu- 
clei. The value C£ s B = -0.050(5) MeV reported in liter- 
ature [22] , is fitted in order to reproduce the difference in 
A-separation energy of the A = 4 mirror hypernuclei ( A H 
and %Re). According to Eq. (6), the contribution of the 
charge symmetry breaking term depends on the differ- 
ence between the number of neutrons and protons. For 
N > Z the CSB term is strictly positive. This implies 
a repulsive contribution per neutron excess that would 
further lower the -Ba for where there are 10 more 

neutrons than protons, thereby reducing the discrepancy 
with the experimental result. 
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Figure 2. A-separation energy as a function of the baryon 
number A. Plain green dots [dashed curve] are the avail- 
able Ba experimental values. Empty red dots [upper banded 
curve] refer to the AFDMC results for the nuclear AV4' po- 
tential plus the two-body AN interaction alone. Empty blue 
diamonds [lower banded curve] are the results with the inclu- 
sion of the three-body hyperon-nucleon force. 

In this letter we have presented the first accurate cal- 
culation of the A-separation energy for closed-shell A- 
hypernuclei using the available microscopic interactions. 
Using the AFDMC method we were able to extend the 
calculation in the medium-heavy range of hypernuclei up 
to A=91, providing, for the first time, a consistent calcu- 
lation of light and heavy hypernuclei. The main outcome 
of the study is that the inclusion of the three-body ANN 
interaction is fundamental in order to reproduce the satu- 
ration properties of the A binding energy in hypernuclei. 
The leading contribution to the three-body interaction 
is strictly repulsive in the range of hypernuclei studied. 
Within the model that we have studied, the inclusion 
of the AN force without a three-body force gives a very 
unphysical A binding energy. 
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We speculate that this would lead to a stiffer EOS for 
the A-neutron matter when the presented interaction is 
applied to the study of the homogeneous medium. This 
fact might eventually reconcile the onset of hyperons in 
the inner core of a neutron star with the observed masses 
of order 2 M . A study along this direction is in progress. 
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